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Abstract. In a recent paper [I], Baxter and Zeilberger show that the two 
most important Mahonian statistics, the inversion number and the major in- 
dex, are asymptotically independently normally distributed on permutations. 
In another recent paper [2], Canfield, Janson and Zeilberger prove the result, 
already known to statisticians, that the Mahonian distribution is asymptoti- 
cally normal on words. This leaves one question unanswered: What, asymp- 
totically, is the joint distribution of the inversion number and the major index 
on words? We answer this question by establishing convergence to a bivariate 
normal distribution. 



1. Introduction 

Let S n be the set of permutations of n objects, that is the set of bijections from 
[n] = {1,2, ... ,n} to itself. Define the inversion number as invfjr) — £ 
[n] x [n] | i < j and ir(i) > 7r(j)} and the major index as J2i£Des(-K) h wnere 
Des(ir) — {i G [n — 1] | > ir(i + 1)} is the descent set of ir. It is a classical 
result due to MacMahon [8] that these have the same distribution, termed the 
Mahonian distribution in his honour. Their probabilty generating function with 
respect to the uniform probability on S n is 

V g to °W p r ( 7r ) = y Pr ( 7r ) = 1 ff 

^— i n\ J -- L l — o 



TreS n 7rGS„ 

There is a bijection <j> from S n to itself due to Foata [1] that satisfies inv((f)(Tr)) = 
maj(n), again proving that inv and maj are equidistributed on S n . This bijection 
fixes the last letter of a permutation, so inv and maj are even equidistributed on 
S n .i — {tt S 5* n | 7r(n) — i} for all i £ [n]. 

In a recent paper [T] , Baxter and Zeilberger show that inv and maj are asymptot- 
ically jointly independently normally distributed on S n , that is 



Pr( m ^<,and ^tUi < y ) J_ / e ^ I ^ 

as n — > oo, 

where u — u n — E(m») = E(maj) is the mean and a 1 = a\ = E((inw — u) 2 ) — 
E((maj — u) 2 ) is the variance of the Mahonian distribution. 

A simple generalization of permutations are words. A word of length n is just 
a function w : [n] — > [d\. It can be regarded as a permutation of a multiset, some- 
times called a permatution. Now the inversion number and the major index can 
be defined for words in the same way as for permutations. Foata's bijection 4> also 
extends to words 0J, so if A is a multiset, and Sa the set of permatutions of it, 
inv and maj are equidistributed on Sa- In fact, as for permutations, they are even 
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equidistributed on S A ,i = {w E S A \ «>(Ej=i a j) = 0- If A = {! Ql ,2 a2 , . . . ,d a "} 
is the multiset containing a\ occurences of 1, a 2 occurences of 2, . . . , a<j occurences 
of d, their probability generating function on Sa is 

fai-\ ha d \ 

J2 q mv{w) Pr(w) = J2 r ajiw) r*M= ):xz!y 

w£Sa w£Sa ^ ai,....a d I 

where f aiH ha<i ) = fer-i "t° d K ! is the g-multinomial coefficient, defined in terms 

V ai,...,a d ) q lai\gl---[adlgi * ' 

of the g-factorials [n] q \ = [l] q ■ ■ ■ [n] q and the g-numbers [n] q = (1 — g")/(l — g). 

In a recent paper [2], Canfield, Janson and Zeilberger use methods from experi- 
mental mathematics as well as Fourier analytic techniques to give multiple proofs 
of the fact, already known to statisticians, that the Mahonian distribution is asymp- 
totically normal on Sa as well. This leaves one question unanswered: What, asymp- 
totically, is the joint probability distribution of inv and maj on Sa? 



2. The Main Result 

Theorem 1. The inversion number and the major index are asymptotically jointly 
normally distributed on words. That is, if A is the multiset containing m\a oc- 
curences of 1, m 2 a occurences of 2, . . . , m^a occurences of d, then ( mv ~f J - ; ma ^~ A1 ) 
on Sa tends to a bivariate normal distribution as a — >• oo. The correlation coeffi- 
cient of that distribution is 

£i<t<j<d "^to? - Ei<i<j<d mjrrij 

Ei<i<j<dOi TO i + m i m j) + 2 T,i<i<j<k<d mm^nik ' 

3. Proof 

To prove this, we use the method of moments. Recall that two random variables 
X and Y jointly converge to a bivariate normal distribution if and only if their 
mixed moments E(X r Y s ) converge to the mixed moments of the bivariate normal 
distribution [3J. We shall instead consider factorial moments M(X-Y-), with the 
factorial powers x— — x(x — 1) • • • (x — k + 1), derive a recurrence for them, and 
check that up to the leading terms it agrees with a well-known recurrence for the 
mixed moments of a bivariate normal distribution, thereby proving the result. 

Consider the effect that removing the last element from a word has on its inversion 
number and its major index. The inversion number will decrease by the number of 
letters smaller than the last letter in the word, and the major index will decrease 
only if the second to last letter was bigger than the last, by the position of the 
second to last letter. In the language of generating functions this is 
(1) 

F(a,i)(p,q) =p^. + i«#F(a-e 1)J )(p,,) + ^=i«n £ F(a- ei , j)(p, q)\ 

\ ]= i 3 =i+i ' 

where F( a ,i)(p,q) = F((a u a 2 , ■ ■ ■ , a d ), i)(p, q) = Y. w eS A , i f nv{w) 1 maj{w) k ^ 
double generating function of inv and maj on the set of permatutions of the multiset 
A = {l ai , 2 a2 , . . . ,d°" d } that end with the letter i, and e* is the i-th unit vector, 
having a 1 in its i-th coordinate and zeroes in all others. The mean of inv on Sa 
is easily seen to be e2(a)/2, where e2(a) = ^ Ki<3<d a^aj is the second elementary 
symmetric polynomial in a-i, . . . ,a,d, and maj has the same distribution, so the 
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centralized probability generating function of inv and maj on Saa is 

F(a,i)(p, q) 



G(a,i)(p,q) 



e 2 ( a ~ e i) I 



where (^i= iai ) = (^ J=1 ^) is the multinomial coefficient. So the recurrence (f) 
translates to 

E ^ °') (pq)^^U + i ^ G (a + ei ,i)(p, q) 



+ '- 1 ^ f^ 1 °*- 1> )(M)^ +E ^-»G(«,i)( P , 9 )Y 

,-=i+i V a e J / / 

which simplifies to 
(2) G(a + e,,z)(p, g ) 

z 

j=i Z)*;=i °fe 

d 

+ ^ pKELj+i "t-Ew «fc) 9 2(^fc= 3+ i »t+EUi o*-E*i+i afc) G(a, q). 

Let FM(&, i, r, s) — E((int)- fi)-{maj — /i)-) be the (r, s)-th mixed factorial moment 
of the random variables inv — fi and maj — n, where \x = E(ira) = E(maj), on 
Sa,». Then 

So we get the Taylor series 



G(a,0(l +P ,l + g )=E fM( Hj' r,a V g '. 

r,s— 

Thus (2) translates to 

E FM(a + ej,i,r,s) r s 
r!s! P 9 

r,s— 
z 

= ^ ^2 (i _|_ p}i(T,i= j+ i a *-Ei=i a k)(i _l_ g)i(E]U J+ i a fc-Ei=i «fc-Efc=i+i a fc) 



■DC 



FM(a,j,r,a) 
r!s! '' ! 



• E 

r,s— 

4- — ^ (i 4- 2 (Efc=j+i a fc~Ei = i a *.0 ^ _|_ ^ 2 (Efc=j+i 

j=i+l Sfc=i a *= 



r!s! 

r,s— 

Now expand the powers using the binomial theorem and compare coefficients to 
obtain 
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FM(a + ej, i, r, s) 

E a i y^ y^ / 2^fc=j+i °fe _ z^fe=i a * 

'|(Efc=i+l a k - Efc=l °fe ^ Efc=i+1 °fc 



r — r' 



y^ fl j y^ y^ 

j=i+l Z^fe=l a fe r'=0 s'=0 

KEfe=j+i °fc + El=i °fe - Efe=l+i 



tL s !=i-FM(a,i,r' )S ') 



- S ^M(a,j,r' (S '). 



So 

d d 
( 3 ) (5^o J -)jJ , Af(a+e i ,t,r,a)-5^a i iJ , M(a,j,r,s) 



l(Efc=,-+i - Efe=i a fc. 



S> E ( 2V ^ + r _ r , 

J=l {(r', s ; ) |r'<r or s'<s} 



-^M(a,i,r',s') 



^ a . /|(Efe=j+i a fe ~ Efe=i a k 

j=i+l {(r' ,s')\r' <r or s'<s} 
Yi(Efc=j+l + Efe=l - Efe=i+1 "fe)\ r -r' s- S ' FA j 



FM(a,j,r', S ')- 

s — s' 

To really profit from this recurrence, we need the following lemma. 

Lemma 1. FM(a,i,r, s) is a rational function in ai,...,a,d, with denominator 

( Ej=i oj - 1) ■ 

Proof. We use a combinatorial argument similar to one given by Baxter and Zeil- 
berger for permutations [T]. For simplicity, we shall first consider ordinary mixed 
moments instead of factorial moments. Start with the equation 

¥ l (inv r maj s ) — 

Ed 
3 = 1 a . 

< k and w(j) > w(k)} 



/ 2^3=1 a-j 

U^~\ E E-' 

a _ ei ) wGS A ,i V j,k=l 

[w(j) > w(j + 1) and k<j] \ , 
j,k=i J 



where [P] = 1 if P is true and [P] — otherwise. Now expand the powers and 
move the sum over Sas and the division by ^j=^°4 -1 ) to the right of all the re- 
sulting 2(r + s) summation symbols. Then there are only finitely many possibilities, 
independent of a±, . . . ,ad, of choosing the relative order of the 2(r + s) indices of 
summation. Some of these may coincide, so say there are I distinct indices i\, . . . , i\. 
Let O = {j 6 [I] | ij + 1 = ij+i}, so o — \0\ is the number of distinct indices whose 
successor is also an index. Now there are only finitely many possibilities, indepen- 
dent of ai, . . . , ad, of choosing the relative order of the indices, then picking the set 
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0, then specifying whether ii = Y,j=i a h H = Sj=i a j ~ 1 or *i < Sj=i a i - 1) 
and afterwards choosing the values taken by w on {ii, . . . , ii} U {i\ + 1, . . . , ii + 1}. 
For each of these scenarios, the summand takes the same value, either or 1, on 
all of the words that belong to that scenario. The number of words belonging to a 

scenario is fef') F U T^) * U = <h> PT-Ti^) ^HX^) 

if i} = Y,U - 1, an d f ^^r^ 1 ) FHX-T' 1 ) if «' < Ej=i a l ~ L where 
v is the vector of multiplicities of the values taken by w on the positions of the 
indices and their successors as defined by the scenario. Here the binomial coeffi- 
cient counts the number of ways to choose the positions of the indices, while the 
multinomial coefficient counts the number of ways of choosing the values of w at 
positions smaller than Ylj—i a j that are not an index or the successor of an index. 
So the contribution to the sum of that scenario is either zero, 



aj-h a 3 -2l+o+U ,T, d 3 = 1 aj-l-U ( ^ =1 a 3 -2i+o\ 
V l-o-l A a-v / V l-o-l A a-v / 



a— e. 



V l — o J V a— v— ej / 

or 



V a— / 



21- 



But all of these are rational functions with denominator {Yli=i a j — -0 > so 



E(inv r maj s ) is a rational function with denominator (X)j=i a j — l) 4 ^ r ~*~ s \ From 
this we can recover E((inv — fj,) r (maj — /i) s ) using the binomial theorem and 
FM(a, i, r, s) = E((inv— /i)— (maj— fj,)-) using Stirling numbers [BJ. So FM(&, i, r, s) 
is a rational function with denominator (^j=i a j — l) 4 r+ . □ 

In fact, computational experiments suggest that the FM(&, i, r, s) are even poly- 
nomials in ai, . . . , ad, but the author has been unable to prove this. 

So now wc know that the FM(a, i, r, s) are of the form d P ^ a ' t,r ' s j (r+a) , where 

P(a, i, r, s) is a polynomial in the variables ax, ... } ad- But what is the degree of 
P(a,i,r,s)? 



Lemma 2. The degree of P(a,i,r, s) is at most 4p(r + s). 

Proof. Let a = (ai,...,ad) be a vector of nonnegative integers and let A = 
, . . . , d ad }. Define X & i — inv — ijl and Y a i = moj — \i on SU.,-*- Then 

|E(^a,)l < E(\X a ,\ r \Y a ^) < E((|X a ,D^)*E((|y a!i r)^)* 

= E(|x aji r+ s ) = e(x;+ s ) 

if r + s is even, by Holder's inequality. From the the paper by Canfield, Janson and 
Zeilberger 0, we know that E(X^+ S ) is a polynomial of degree |(r + s), so that if 

E(X^ )i l^ )i ) is written as a rational function with denominator ( Ylj=i a j ~ l) 
the degree of the numerator is at most 4(r + s) + 1 (r + s) = ^-(r + s). FM(a,i,r, s) 
can be recovered from the E(X^ i Y a s i ) using Stirling numbers [BJ, so the degree of 
P(a, i, r, s) is also at most 4r (r + s). If r + s is odd, 

< E(|X a)i r+ s ) < E(\X & , i \ r + s + 1 )^ = E(X r a +° +1 )^h 
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by Jensen's inequality, so the degree of P(a, i, r, s) is at most 

4 ( r + *) + f( r + 8 + 1 ) r +t+i = T (r + s) - 



□ 



Remark If r + s is odd, ^(r + s) is not an integer, so the degree of P(a, i, r, s) 
is strictly less than ^-(r + s). 

For r+s even, we substitute FM(a, i, r, s) = P< ' a ''' r \ s i (r+3 ) in (3) and multiply 
by the common denominator ( X^jLi a j — l) 4 ^'^ to obtain the polynomial identity 



a a 

(4) -4(r + sj^P(a + e i ,i,r,s)-^2a j P(a,j,r,s) 



3 = 1 



d 



5^ afe ) 2! 51 afe ~ 5Z afc ) rP ( a '- ? ' r ~ 1 ' s ) 

fe=i fe=j+i fc=i 

4 2 * 17—1 d 

+ 2! 51 a k -J2a k - a k ^sP(a,j,r,s - 1) 

fc=l fe=j+l fc=l fe=i+l 

+ (5> fc ) 2 ( 2 ( 5^ afc -E afc )j Kr-l)P(a,j,r-2,s) 
fc=i ^ k=j+i fe=i ' 

d g ^ i j — 1 d s 2 

+ (5Z flfc ) 2 ( 2( 5^ a fe"5Z afc_ 5Z afc )j s(s ~ !) p (a,i>^s - 2 ) 

fe=l ^ k=j+l fe=l fc=»+l ' 

d g j d j — 1 i j — 1 d 

+ (5I afe ) 2! 51 ^"X^)^ 51 a k-^2a k - Y a k ^rsP(a,j,r-l,s-l) 



k=j+i k=i 



+ a J 

j=i+l 
d 



k=j+l fe=l fe=i+l 



(5Z afc ) 2! 51 afe ~ 5I a ' £ ) rjP ( a '- ? ' r ~ 1 ' s ) 



fc=i+i 



fc=i 
i-i 



+ 



fe=i 



fc=i+i 



fe=l fe=i+l 
d 



+ (5I afe ) 2! 51 o-k + ^a-k- a k ^sP(a,j,r,s-l) 

k=l k=j+l fc=l fe=i+l 

d 8 ]/ 1 d J_ 1 \ 2 

(5Z° fc ) 2 ( 2 ( 5^ flfc ~5I afc )) ^(^-l)-P(a,j,r-2,s) 
fe=i ^ fe=j+i fe=i ' 

d 8 1 / 1 ^ ' ^ 1 \ 2 

+ (5I afe ) 2 ( 2 ( 5^ a k + Y ak ~ 51 afc )j s ( s - l ) p (^^ r ^ s - 2 ) 



+ (5I afe ) 2! 51 ^"X^)^ 51 a fc+X] afe ~ 51 a fc)rsP(a,j,r-l,s-l) 

fe=l k=j + l k=l k=j+l k=l k=i+l 

+ lower order terms. 
Here the lower order terms have degree at most ^(r + s) — 1. 
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Let us now set up a recurrence that is satisfied by the mixed moments of the 
bivariate normal distribution (X,Y), with E(X) = E(Y) = 0, 

Var(X) = Var(Y) = (a;a 2 + afaj) + 2 didjCtk and 

l<i<j<d l<i<j<k<d 

Cov(X,Y) = ctiOj — ajcij. 

l<i<j<d l<i<j<d 

Recall that the mixed moments of a multivariate normal distribution are given by 
Isserlis' Theorem [7J: 

E(X 1 ---X„) = ^ JJ EpQX,), 

P 

where the sum is over all perfect matchings of [n] , which holds if X± , . . . , X n are 
jointly normally distributed with zero mean. In this case, this means that 

(5) E{X r Y s ) = ^E(X 2 )l p lMlE(y 2 )l p lwlE(Zy)l p|K ".«l, 

p 

where the sum is over all perfect matchings of [r] II [s], P|[ r ] is the set of edges of 
P within [r], P|[ s ] is the set of edges of P within [a], and P|ic is the set of edges 
of P that have a vertex in [r] and [s] each. Now double count perfect matchings 
of [r] II [s] with the above weighting and a single vertex coloured red: On the one 
hand, one can pick a perfect matching and then choose a vertex to colour red. On 
the other hand, one can pick a vertex to colour red, then choose a vertex to match 
to it, and afterwards take a perfect matching of the rest. This leads to the formula 

(r + s)E(X r Y s ) = E{X 2 )r(r - l)E{X r - 2 Y s ) 

+ E(Y 2 )s(s - l)E{X r Y s - 2 ) 

+ E(XY)2rsE{X r ~ 1 Y s ~ 1 ). 

Now that we have these two recurrences, we can prove the following fact. 

Claim 1. For r + s even, the leading terms of P(a, j,r, s) are the same as those of 
(Xw=i a j) 4 r+S ^(X r Y s ), where (X,Y) is distributed as above, for all i 6 [d]. 

Proof. We proceed by induction on r + s. The result is true for r = s = 0. If 
r + s > 0, consider the recurrence (4). Multiply this equation by at and sum 
over all i £ [d], simplifying the right-hand side using the fact that, by induction 
hypothesis, the leading terms of P(a, j, r' , s') do not depend on j for r' + s' < r + s. 
After some amount of routine calculation, the result is 

d , d d \ 

(6) X! a M (5Z a J ~ 4(r + s)JP(a + ei,i,r,s) - ^OjP(a, j, r, s) J 

i=l V 3 = 1 3 = 1 ' 

\ d 9 

= g(X! afe ) ( X! (a 2 a 3 +aia 2 ) +2 ^ OtOjafcJ r(r — l)P(a, J,r — 2, s) 

k=l l<i<j<d l<i<j<k<d 

\ d 9 

+ g(X! a ' c ) ( X! (a 2 aj + a,a 2 ) + 2 ^ ctiajak)s(s -l)P(a,l,r,s -2) 



k=l l<i<j<d l<i<j<k<d 

d 



jjr(X) a fc) ( X! a * a J ~ X! OiOj) 2 ^sP(a,Z,r-l,s-l) 

+ lower order terms, 



k=l l<i<j<d l<i<j<d 
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for any I G [d]. Notice that all the terms involving P(a, i,r — l, s) or P(a, i,r,s — 1) 
conveniently cancel out. 

Now, by induction hypothesis, the right-hand side of (6) has degree -^(r + s) + 1. 
Suppose for contradiction that the degree of P(a, i,r, s) is less than ^-(r + s), for 
all i G [d]. The leading terms of the left-hand side cancel out, so the degree of 
the left-hand side is less than -^-(r + s) + 1, which is absurd. Thus there exists 
io G [d] such that the degree of P(a, i , r, s) is exactly -y (r + s). The degree of the 
right-hand side of (4) is at most -y(r + s), thus the leading terms of the left-hand 
side have to cancel out for i = i - By substracting (4) for i = i from the same 
equation for i = i\ ^ i , we see that the leading terms of P(a, i, r, s) do not depend 
on i, so in particular the degree of P(a, i,r, s) is ^-(r + s) for all i G [d]. Recall 
that, for a polynomial /(a) = /(oi, . . . , a^), 

a 

/(a + e^) = /(a) + — — /(a) + lower order terms 

OQ>i 



and 

da 



^ai— — /(a) = deg(/)/(a) + lower order terms. 



From these two observations, and (6), there results 

( E °0 ( y( r + s ) - 4 ( r + *)) p ( a ' r ' s ) 

i— 1 ^ ' 
d 9 

= s(E afe ) ( E K 2 aj+aja|)+2 ^ a^a^ 

fe=l l<i<j<d l<i<j<fe<d 

• sr(r - l)P(a,Z,r - 2,s) 

\ d 9 

+ s(E afe ) ( E (afa 3 +a,a 2 J ) +2 ^ a^a^ 

fc=l l<i<j<d l<i<j<k<d 

■ s(s-l)P(a,l,r,s-2) 

\ d 9 

+ 8(E flfe ) ( E a * a ? ~ E otaj)lrsP{a.,l,r-l,s-l) 

k=l l<i<j<d l<i<j<d 

+lower order terms, 

for all I G [d]. This gives that, for all I G [d], 
(r + s)P(a, Z, r, s) 

fc=l l<i<j<d l<i<j<fe<d 

•r(r-l)P(a,Z,r-2,s) 

fe=l l<i<j<d l<i<j<fe<d 

■s(s- l)P(SL,l,r,S-2) 

-yd g 

+ 12(E afc ) ( E a * a J ~ E a 4 2 aj)2rsP(a,/,r-l,s-l) 

fe=l l<i<j<d l<i<j<d 

+lowcr order terms. 
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By induction hypothesis, this is further equal to 



k=l l<i<j<d l<i<j<k<d 

■r(r-l)£m Epr- 2 Y S ) 
fe=l 

fe=l l<i<j<d l<i<j<fc<d 

/jL \4(r+s-2) 



fe=l l<i<j<d l<i<j<c 

4(r-l+s-l) 



/ \ iir-ii-s-ii 



fc=i 

+lower order terms. 



(r + s) f 2J Ofe) E(JTY S ) + lower order terms. 

fc=i 



□ 



To conclude the proof of the main theorem, if m = (mi, . . . , mj) is a vector of 
nonnegative integers, and A = {l"' 1 ", 2" i2a , . . . , d mda }, then, on S A)i , 



E 



mv 



[i\ r ( maj — n\ S \ E((inv — /j,) r (maj — /i) s ) 

T r+s 



a' 



—J {— 

FM(am,i,r,s) f1 . E(X r Y s ) ,. _ ( ( X\ r (Y\ s \ 

for r + s even and (X, Y) distributed as above.. For r + s odd, the degree of 
P(am, i, r, s) is less than -^-(r + s), so 

inv — ^mcij — fi^ s ^ E((inv — fi) r (maj — fi) s ) 

= + *> - o - *> = - ( (i ) r © i - *> 

Now (f , docs not depend on a, so E((^) r (22pi) a ) e^)^)') f or 
all r, s. Thus (^^, ^^-^) on S A ,i tends to (f , in distribution as a -> oo. 
A fortiori, the same is true on SU- This is a bivariate normal distribution, with 
correlation coefficient 

E (X Y\ _ E{XY) _ Ei<i<j<d m i m j - Ei<i<f™2™< 



crcrj cr 2 Ei< 4<J <d( m * TO | + m i m i) + 2 Ei<i<j<k<d m l m ] m k ' 



This completes the proof. 
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4. Discussion 
Corollary 1. IfCov(X,Y) = 0, that is 

mim 2 — mfnij = 0, 

l<i<j<d l<i<j<d 

then inv and maj are asymptotically independent. This happens in particular if all 
the mi are equal. 

Also consider the special case d — 2. In this case, the leading terms of -P(a, i, r, s) 
have a particularly simple form. To see this, recall from (5) that 

E(X r Y s ) = J2e(X 2 )\ pi m\e(Y 2 )\ pi m\e(XY)\ PIk ^\ 
p 

r 

= J2 E(x 2 )^-m(xYy 

S=0 {r\\P\«r..H} 

Now if d = 2, Var(X) = Var(Y) = ^gh£ an d Cov(X,Y) = so for r + s 

even, this is equal to 

v f (r-i - m . - 1 - m ->V 

{j£[r]\r—jeven} 

r+s 

ab 



r-(r-s-l)\\{-) (b-a) s (a + b)- 



2 



F 1 - s /2,- S /2+l/2;(r-s)/2 + l; 



a + 6 x 2> 



b — a 



,^,-.s-:i,!!(g > ) ( - L >„n„ + /,r 



• 2 F X ( -s, (r - s)/2 + 1/2; r - s + 1; 



r+s 

ab 



= (12) (a+&)^fo s (r + .s-l)!! 



2 Fi(-s, (r - s)/2 + 1/2; -(r + s)/2 + 1/2; -a/6) 



ab 



I2) 



„ . \ I I • '■ ^~ I .1 1 • • \ 1 ■ '■ ^~ 1 .1. J • - \ !(■//' 

(V - ,s - 1M! v J ' v J ; v ; 

Here we used the hypergeometric identities (JS], (5.10)) ([ID], (1.8.10)) 



2 Fi(a,a + l/2;6+l/2;z 2 ) = (1 - z)- 2a 2 F x ^2a, 6; 26; "-^yj , and 
2-^1(0, -n; c; 2) = (1 — z)™-^^- 2 -F 1 i(c - a, -n; 1 - a - n; 1 - z), 

(Cjn 

and assumed, without loss of generality, that r > s. In particular, the degree of 
¥,(X r Y s ) is |(r + s), a fact already implicitly used in the proof of Claim 1. 
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